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摘 要 : 本 文 提出 了 一 种 通过 在 叶片 前 缘 前 设置 微小 圆柱 来 抑制 叶片 上 流动 分 离 的 控制 方法 ， 并 对 采用 此 流动 控制 方法 后 的 升 
力 型 垂直 轴 风 力 机 的 气动 性 能 进行 了 数值 研究 。 结 果 表 明 ， 在 大 攻 角 情况 下 ， 辟 型 前 缘 前 设置 微小 圆柱 可 以 有 效 的 抑制 风力 机 
叶片 上 的 流动 分 离 ， 得 到 增 升 减 阻 的 效果 同时， 与 传统 的 升力 型 垂直 轴 风 力 机 相 比 ， 前 缘 前 设置 微小 圆柱 后 的 垂直 轴 风 力 机 
在 低 尖 速 比 下 的 风能 利用 率 可 得 到 大 幅度 的 提高 。 
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Abstract: This paper presents a flow control method using micro-cylinder installed in front ofthe blade leading edge 
and invests the aerodynamic performance of the VAWT equipped with micro-cylinders. Numerical results suggest 
that the flow separation ls effectively controlled by setting a micro-cylinder in front of the blade leading edge when 
the blade works at large attack angle. Meanwhile, the VAWT equipped with micro-cylinders has higher efficiency 
than regular VAWT under small tip speed ratio. 
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efficiency 
前半 设置 微小 圆柱 的 流动 控制 方法 ， 并 通过 课题 组 自主 
村 研发 的 计算 软件 UCFDDI 对 采用 这 种 方法 来 抑制 重 


风能 是 一 种 无 污染 的 可 再 生 能 源 ， 它 是 可 再 生 直 轴 风力 机 叶片 上 流动 分 离 的 效果 进行 了 数值 模 
能 源 中 发 展 最 快 的 清洁 能 源 ， 同 时 风力 发 电 也 是 最 拟 。 同 时 ， 对 在 叶片 前 缘 前 设置 微小 圆柱 后 的 垂直 
有 大 规模 开发 和 商业 化 发 展 前 景 的 发 电 方式 。 如 轴 风 力 机 的 气动 性 能 进行 了 数值 模拟 计算 ， 并 将 得 
果 风 能 的 1% 被 利用 ， 则 可 以 减少 世界 3% 的 能 源 消 到 的 风能 利用 率 与 传统 的 垂直 轴 风 力 机 的 风能 利用 
耗 ; 风能 用 于 发 电 , 可 产生 世界 总 电量 的 8%~9%D]。 率 进行 了 对 比 。 
质量 轻 、 结 构 简 单 、 平 衡 性 好 、 安 装 维护 方便 的 s a 
Darrieus 直 叶 片 垂直 轴 风 力 机 (SB-VAWT) 具有 1/ 1 计算 模型 
阔 的 应 用 前 景 所 。 但 传统 的 升力 型 垂直 轴 风 力 机 在 计算 中 升力 型 垂直 轴 风 力 机 叶片 翼 型 选用 
低 尖 速 比 下 常 处 于 大 攻 角 工 况 下 运行 ， 随 之 产生 大 NACA0015, 翼 型 弦 长 c=400 mm, 旋转 直径 D=2500 
尺度 的 流动 分 离 ， 导 致 翼 型 的 升力 系数 降低 ， 阻 力 mm， 自 由 来 流速 度 二 10 m/s 由 。 图 1 为 该 风力 机 的 
系数 增 大 ， 最 终 使 风力 机 叶片 的 气动 性 能 恶化 。 轴 向 截面 示意 图 ， 其 中 0 为 旋转 中 心 ，R 为 旋转 半 

为 了 能 使 升力 型 垂直 轴 风 力 机 在 低 尖 速 比 下 的 径 。 图 2 为 叶片 前 缘 前 设置 微小 圆柱 后 的 放大 图 ， 
气动 性 能 有 所 提高 ， 本 文 提 出 了 一 种 在 叶片 前 缘 前 
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工 为 前 缘 前 微小 圆柱 距离 叶片 前 缘 的 距离 ，d 为 前 
缘 前 微小 圆柱 的 直径 。 


图 1 垂直 轴 风 力 机 的 轴 向 截面 示意 图 
Fig.1 Schematic of the cross section of VAWT 


图 2 叶片 前 缘 前 设置 微小 圆柱 示意 图 


Fig.2 Schematic of the micro-cylinder in front of blade leading 


edge 
程序 基于 三 维 结构 化 网 格 ， 整 个 计算 域 采用 结 
构 网 格 ， 外 边界 距 咽 型 弦 长 中 心 为 80 c， 尽 量 减 小 
边界 对 风力 机 流 场 的 影响 。 同 时 ， 对 叶片 和 微小 圆 
柱 近 壁 面 进行 了 加 密 处 理 。 图 3 为 整个 计算 域 网 格 ; 
4 为 器 型 附近 网 格 。 
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图 3 整个 计算 域 网 格 
Fig.3 Whole mesh of VAWT 


图 4 翼 型 附近 网 格 
Fig.4 Mesh around the airfoil 

2 数值 方法 

本 文通 过 求解 二 维 不 可 压 非 定常 N-S 方程 来 数 
值 模拟 风力 机 流 场 ， 选 用 Spalart-Allmaras (S-A) 注 
流 模 型 ， 空 间 采 用 Roe’s 三 阶 迎 风格 式 ， 计 算 域外 
边界 为 圆 形 ， 边 界 条 件 设 置 为 速度 进出 口 ， 姻 型 和 
微小 圆柱 的 边界 条 件 设 为 无 滑 移 壁面 ， 流 场 来 流 x 
方向 马赫 数 M。 为 0.0294。 


3 计算 结果 及 分 析 


3.1 前 缘 前 设置 与 未 设置 微小 圆柱 时 翼 型 的 气动 性 
能 的 比较 
3.1.1 模型 验证 

为 了 验证 数值 模拟 的 可 靠 性 ， 本 文 首先 模拟 了 
单个 NACA0012 辟 型 在 雷诺 数 Re 为 6000000 的 情 
况 下 ， 该 经 型 升 阻力 系数 随 来 流 攻 角 改变 而 变化 的 
情况 ， 并 将 计算 结果 与 修正 后 的 实验 值 忆 进行 了 比 
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较 。 图 5 为 不 同 攻 角 下 ， 辟 型 升 阻力 系数 的 计算 值 
与 实验 值 的 对 比 图 。 


Coefficient 


angle of attack 


5 不 同 攻 角 下 升 阻 力 系数 对 比 图 
Fig.5 Comparison of the lift and drag coefficient between 


experiment and numerical Simulation at various attack angle 

从 图 5 中 可 以 看 出 ， 攻 角 在 17 "之 前 升 阻力 系 
数 与 实验 值 都 拟 合 得 较 好 ， 但 在 大 攻 角 下 出 现 了 一 
定 的 偏差 。 这 与 数值 模拟 时 采用 二 维 非 定常 模拟 ， 
忽略 了 实际 的 三 维 效应 有 关 ， 且 在 大 攻 角 情况 下 ， 
叶片 上 出 现 了 流动 分 离 ， 流 动 变 得 复杂 。 但 从 总 体 
上 来 看 ， 本 文 所 采用 的 网 格 与 计算 程序 用 于 该 类 咽 
型 的 模拟 还 是 可 靠 的 。 
3.1.2 前 缘 前 设置 与 未 设置 微小 圆柱 时 收 型 的 气动 
性 能 的 比较 

影响 采用 微小 圆柱 抑制 囊 型 上 流动 分 离 效果 的 
主要 因素 有 : 微小 圆柱 的 直径 d 和 微小 圆柱 距离 必 
型 前 缘 距 离 工 。 因 此 ， 本 文 首 先 对 微小 圆柱 的 直径 
ad 和 微小 圆柱 距离 辟 型 前 缘 距 离 工 进 行 了 初步 优化 。 
当 来 流速 度 为 10 m/s， 攻 和 角 18 "时 得 到 的 数值 模拟 
结果 如 图 6 所 示 。 
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图 6 不 同 d 和 工 的 升 阻 比 数值 模拟 结果 
Fig.6 Lift-drag ratio at various L and d 
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从 图 6 可 以 看 出 ， 在 攻 角 为 18* 的 工 况 下 ， 微 
小 圆柱 距离 叶片 前 缘 距 离 工 为 0.03 m 和 0.04m 时 ， 
微小 圆柱 直径 的 增 大 ， 可 使 经 型 的 升 阻 比 降低 ， 同 
时 微小 圆柱 直径 d 为 0.02 m 和 0.03 m 时 ， 上 距离 翼 
型 前 缘 太 近 或 太 远 都 可 导致 翼 型 的 升 阻 比 降 低 。 适 
当 的 微小 圆柱 的 d 和 工 能 使 辟 型 的 升 阻 比 提 高 ， 不 
适当 的 微小 圆柱 的 d 和 工 使 翼 型 的 升 阻 比 下 降 。 表 
1 给 出 了 在 来 流风 速 为 10 m/s， 攻 角 为 18 "情况 下 ， 
前 缘 前 未 加 微小 圆柱 时 的 NACA0012 翼 型 的 升 阻 
力 系 数 与 升 阻 比 , 并 与 在 该 翼 型 前 加 直径 为 0.01 m， 
距 辟 型 前 缘 0.03 m 的 微小 圆柱 后 辟 型 的 升 阻 力 系 
数 和 升 阻 比 进行 对 比 。 

表 1 设置 微小 圆柱 前 后 的 升 阻力 系数 和 升 阻 比 对 比 
Table 1 Comparison of Cl Ca and Ci/Ca between NACA0012 

and NACA0012 with micro-cylinder 


CI Ca CCa 
NACAO0012 0.8928 0.1500 5.95 
设置 微小 圆柱 1.4481 0.09338 15.51 


从 表 1 可 以 看 出 : 前 缘 前 加 了 微小 圆柱 后 的 经 
型 ， 其 升力 系数 可 提高 60% 以 上 ， 且 阻力 系数 有 所 
降低 ， 使 得 翼 型 的 升 阻 比 得 到 提高 。 图 7 和 图 8 为 
在 18 " 攻 角 下 , 在 NACA0012 辟 型 前 缘 前 适当 位 置 
处 未 设置 和 设置 了 微小 圆柱 后 惨 型 周围 流 线 变化 的 
比较 。 


Se Si 3 
图 7 咽 型 前 缘 前 未 设置 微小 圆柱 时 的 流 线 图 
Fig.7 Streamlines of NACA0012 at 18 ° attack angle 


8 辟 型 前 缘 前 设置 微小 圆柱 时 的 流 线 图 
Fig.8 Streamlines of NACA0012 with micro-cylinder at 18 


丽 


attack angle 
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可 以 看 出 ， 在 器 型 前 缘 前 设置 微小 圆柱 后 ， 跟 
型 上 原 有 的 较 大 的 流动 分 离 泡 (如 图 7 所 示 ) 变 成 
了 靠近 辟 型 尾 缘 处 的 很 小 的 分 离 泡 〈 如 图 8 所 示 )， 
轻型 上 的 流动 分 离 得 到 了 显著 抑制 。 

在 以 上 数值 模拟 结果 的 基础 上 ， 试 图 将 这 种 流 
动 控制 的 方法 进一步 应 用 到 升力 型 垂直 轴 风 力 机 的 
叶片 上 ， 当 风力 机 叶片 处 于 大 攻 角 工 况 运 行 时 ， 采 
用 这 种 方法 可 以 有 效 推迟 或 抑制 叶片 上 的 流动 分 
离 ， 提 高 风力 机 叶片 的 气动 性 能 。 

3.2 风力 机 叶片 前 缘 前 设置 微小 圆柱 后 的 气动 性 能 
分 析 及 与 传统 风力 机 结果 的 比较 

在 模拟 叶片 前 缘 前 设置 微小 圆柱 的 升力 型 垂直 
轴 风 力 机 之 前 ， 本 文 先 对 文献 4 中 的 升力 型 垂直 轴 
风力 机 进行 数值 模拟 。 模 拟 结果 与 实验 值 相 比 ， 最 
佳 风 能 利用 率 所 对 应 的 叶 尖 速 比 基 本 一 致 ， 但 最 佳 
风能 利用 率 偏 高 。 风 能 利用 率 偏 高 与 数值 模拟 时 采 
用 二 维 非 定常 模拟 ， 忽 略 了 实际 三 维 效应 有 关 。 综 
上 ,本文 所 采用 的 网 格 与 计算 程序 是 有 效 的 。 接着 ， 
本 文 对 叶片 前 缘 设置 微小 圆柱 的 风力 机 的 气动 性 能 
进行 数值 模拟 。 

影响 采用 微小 圆柱 抑制 风力 机 叶片 上 流动 分 离 
效果 的 主要 因素 有 : 微小 圆柱 的 直径 d 和 微小 圆柱 
距离 叶片 前 缘 距 离 工 。 因 此 ， 本 文 对 微小 圆柱 的 直 
径 d 和 微小 圆柱 距离 风力 机 叶片 前 缘 距 离 工 进 行 了 
初步 优化 。 当 来 流速 度 为 10 m/s， 叶 尖 速 比 为 1.1 
时 ， 采 用 不 同 d 和 工时 的 风力 机 风能 利用 率 的 模拟 
结果 如 图 9 所 示 。 
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9 叶 尖 速 比 1.1 时 不 同 设置 微小 圆柱 工 况 的 C 对比 图 
Fig.9 Wind energy utilization efficiency of various L and d at 
1.1 tip speed ratio 

从 图 9 可 以 看 出 , 随 着 微小 圆柱 直径 d 的 增加 ， 
风力 机 的 风能 利用 率 有 所 下 降 ， 但 随 着 距离 工 的 增 
加 ， 风 能 利用 率 则 在 不 断 增 加 。 图 10 给 出 了 距离 工 
为 0.02 m， 直 径 d 为 0.001 m 的 微小 圆柱 设置 与 未 
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设置 微小 圆柱 的 风力 机 Cs 在 低 叶 尖 速 比 下 的 变化 


| 线 。 
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图 10 设置 小 圆柱 前 后 Cr 随 叶 尖 速 比 变化 图 


Fig.10 Comparison of wind energy utilization efficiency 


between regular VAWT and VAWT with micro-cylinder at 


various tip Speed ratio 


从 图 10 中 可 以 看 出 ,在 低 尖 速 比 1.1 的 情况 下 ， 


叶片 前 缘 前 设置 微小 圆柱 后 的 风力 机 Cs 比 传统 风 
力 机 提高 了 9%。 在 低 尖 速 比 情况 下 ， 前 缘 前 设置 
微小 圆柱 对 风力 机 气动 性 能 的 提高 是 非常 显著 的 。 
图 11 和 图 12 分 别 给 出 了 在 叶 尖 速 比 1.1 时 风力 机 
叶片 前 缘 前 适当 位 置 处 未 设置 和 设置 了 微小 圆柱 的 


流 线 图 。 


i 
be x 


i s 
图 11 某 时 刻 风 力 机 叶片 前 缘 前 未 设置 微小 圆柱 流 线 图 
Fig.11 Streamlines of regular VAWT 
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图 12 菜 时 刻 风力 机 叶片 前 缘 前 设置 微小 圆柱 流 线 

Fig.12 Streamlines of VAWT with micro-cylinder 
从 图 11、12 中 可 以 看 出 , 在 左下 角 位 置 , 通过 
在 叶片 前 缘 前 设置 微小 圆柱 ， 叶 片上 原 有 的 大 流动 
分 离 泡 被 抑制 , 从 而 改善 了 风力 机 叶片 的 气动 性 能 ， 
提高 了 升力 型 垂直 轴 风 力 机 在 低 尖 速 比 下 的 风能 利 
用 率 。 


本 文 针 对 风力 机 叶片 在 大 攻 角 下 运行 易 产 生 流 
动 分 离 的 特征 ， 提 出 了 一 种 在 风力 机 叶片 前 缘 前 设 
置 微小 圆柱 来 提高 风力 机 叶片 气动 性 能 的 流动 控 珊 
方法 ,并 采用 课题 组 自主 研发 的 计算 软件 UCFD 对 
该 流动 控制 方法 的 有 效 性 进行 了 数值 模拟 研究 ， 得 
到 以 下 结论 : 

(1) 在 大 攻 角 情况 下 ， 在 如 型 前 缘 前 设置 微小 圆柱 
可 以 有 效 的 抑制 辟 型 上 的 流动 分 离 ， 起 到 增 升 减 阻 
的 作用 ; 


EE 
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(2) 与 传统 垂直 轴 升 力 型 风力 机 相 比 ， 在 低 尖 速 比 
下 ， 叶 片 前 缘 前 设置 微小 圆柱 后 的 升力 型 垂直 轴 风 
力 机 的 风能 利用 率 有 明显 提高 , 如 在 叶 尖 速 比 为 1.1 
的 工 况 下 ,在 风力 机 叶片 前 加 设 优 化 后 的 微小 圆柱 
风力 机 的 风能 利用 率 可 提高 9%。 
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